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Abstract: We compute the decoupling relations for the strong coupling, the light quark 
masses, the gauge-fixing parameter, and the light fields in QCD with heavy charm and 
bottom quarks to three- loop accuracy taking into account the exact dependence on mc/rrH). 
The application of a low-energy theorem allows the extraction of the three-loop effective 
Higgs-gluon coupling valid for extensions of the Standard Model with additional heavy 
quarks from the decoupling constant of as ■ 
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1. Introduction 

QCD where all six quark flavours are treated as active degrees of freedom is rarely used in 
practical applications. If the characteristic energy scale is below some heavy-flavour masses, 
it is appropriate to construct a low-energy effective theory without those heavy flavours. 
The Lagrangian of this theory has the same form as the one of QCD plus corrections 
suppressed by powers of heavy-quark masses. Usually, heavy flavours are decoupled one at 
a time which results in a tower of effective theories, each of them differ from the previous 
one by integrating out a single heavy flavour. The parameters of the Lagrangian of such 
an effective low-energy QCD (^^(/x), the gauge flxing parameter a(/x), light-quark masses 
mj(/u)) are related to the parameters of the underlying theory (including the heavy flavour) 
by so-called decoupling relations. The same holds for the light fields (gluon, ghost, light 
quarks) which exist in both theories. QCD decoupling constants are known at two- [Q, ^ |^, 
three- |Q] and even four- loop order Q, g]. 

The conventional approach just described ignores power corrections in ratios of heavy- 
quark masses. Let us, e.g., consider the relation between as and as (the superscript 
denotes the number of active flavours). Starting from three loops, there are diagrams 
containing both b- and c-quark loops which depend on mc/rrib- The power correction 
^ {as/n)^ {mc/mi,)'^ is not taken into account in the standard approach, although, it 
might be comparable with the four-loop corrections of order (os/vr)^. In the present paper, 
we consider {rric/rrn,)^ power corrections at three loops by decoupling b and c quarks in a 
single step. 

Of course, the results presented in this paper are generic and apply to any two flavours 
which are decoupled simultaneously from the QCD Lagrangian. Our full theory is QCD 
with ni light flavours, ric flavours with mass rric, and nb flavours with mass nib (in the real 
world ric = Ub = 1). Furthermore we introduce the total number of quarks nj = ni+ric+nb. 
We study the relation of full QCD to the low-energy effective theory containing neither b 
nor c. 

The bare gluon, ghost and light-quark fields in the effective theory are related to the 
bare fields in the full theory by 

4-) = {c\f"At^-^ , 4-^ = {Cf\t''^ , czS"'^ = {c',f\^:'^ , (1.1) 

where the bare decoupling constants are computed in the full theory via [§ 

C^(«io'^«^^) = 1 + nA(o) = [zir\ 
C°(«lo'\a^^) = i + ne(o) = [zr]-\ 

C°(«l^^\a^^) = l + Sy(0)=[^°r'' (1-2) 

with Q^o = 5o/(47r)^-^; UAiq"^), Udq"^) and S(g) = iT.v{q^) + nigo^siq'^) are the (bare) 
gluon, ghost and light-quark self-energies (we may set all light-quark masses to in Sy and 
S5). The fields renormalized in the on-shell scheme coincide in both theories; therefore. 



the bare decoupling coefficients (L2) are the ratio of the on-shell renormalization constants 



of the fields. In the effective theory all the self-energies vanish at q = (they contain no 



scale), and the on-shell Z factors are exactly 1. In the full theory, only diagrams with at 
least one heavy-quark loop survive.^ 

Next to the fields also the parameters of the full and effective QCD Lagrangian are 
related by decoupling constants 

J^l) _ aO {«/) M _ M(nf) ini) _ ("/) ,. o^> 

where a is the gauge parameter defined through the gluon propagator 

D,Ak) = -^(9,.-{l-a)^) . (1.4) 



^2 y^l^" V" "-^ ^2 
The bare decoupling constants in Eq. ( |1.3| ) are computed with the help of [^ 

C(«io'^) = (1 + ^A-ccf {Zff Z^l = (1 + VAqqf [Z^^f Z^l = (1 + Vaaa? {Z^lf , 
C(ai7^) = Z-[l-S5(0)]. (1.5) 

The Ace, Aqq and AAA proper vertex functions are expanded in their external momenta, 
and only the leading non-vanishing terms are retained. In the low-energy theory they get 
no loop corrections, and are given by the tree- level vertices of dimension-4 operators in the 
Lagrangian. In full QCD (with the heavy flavours) they have just one colour and tensor 
(and Dirac) structure, namely, that of the tree-level vertices (if this were not the case, the 
Lagrangian of the low-energy theory would not have the usual QCD form^). Therefore, we 
have the tree-level vertices times (1 + Fj), where loop corrections Tj contain at least one 
heavy-quark loop. The various versions in the first line of Eq. ( |1.5[ ) are obtained with the 
help of the QCD Ward identities involving three-particle vertices. In our calculation we 
restrict ourselves for convenience to the ghost-gluon vertex. Note that the gauge parameter 
dependence cancels in ^^^ and Cm whereas the individual building blocks in Eq. (ll.5|) still 
depend on a. This serves as a check of our calculation. 

The MS renormalized parameters and fields in the two theories are related by 

c^tHp') = CaAp',fi)a'^'\f^), a("')(/i') = a(/u',;u)a("/)(^), 
4"'^(/^') = U(/x',^)m^)(^) , A("')(/.') = Cy'(/.',/x)A("/)(/.) , 
cK)(^') = cy2(^',^)cK)(^) , g("')(;u') = C,'/2(/i',^)g("^)(/u) , (1.6) 

where we allow for two different renormalization scales in the full and effective theory. 
The finite decoupling constants are obtained by renormalizing the fields and parameters in 



^At low q ^ 0, the self-energies in the full theory are given by sums of contributions from various 
integration regions, see, e.g., B; the contribution we need comes from the completely hard region, where 
all loop momenta are of order of heavy-quark masses. 

^The Aqq vertex at 0-th order in its external momenta obviously has only the tree-level structure. For the 
Ace vertex at the linear order in external momenta, this statement is proven in AppendixM. The AAA vertex 
at the linear order in its external momenta can have, in addition to the tree-level structure, one more struc- 
ture: d''i''2a3(gfiM2;j^3_^cycle); however, the Slavnov-Taylor identity {T{d''A^{x),d''A^{y),d^Axiz)}} = 
leads to T^^imiK"- k!^^ k^'^ = (see Ref. 0), thus excluding this second structure. 



Eqs. (U) and (ITj) which leads to 
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(1.7) 



where Z^ ^ are the MS renormahzation constants in nj-flavour QCD which we need up 
to three-loop order. 



2. Calculation 

Our calculation is automated to a large degree. In a first step we generate all Feynman 
diagrams with QGRAF [Q]. The various diagram topologies are identified and transformed 
to FORM Q with the help of q2e and exp |10, ^ (these topologies have been investigated 
in [^]). Afterwards we use the program FIRE [^] to reduce the two-scale three-loop 
integrals to four master integrals which can be found in analytic form in Ref. |Q]. 

As a cross check we apply the asymptotic expansion (see, e.g., Ref. Q) in the limit 
rric ^ rrif, and evaluate five expansion terms in (mc/mi,)'^. The asymptotic expansion 
is automated in the program exp which provides output that is passed to the package 
MATAD [^] performing the actual calculation. 

In the following we present explicit results for the two-point functions and T^cc needed 
for the construction of the decoupling constants. Other vertex functions can be easily 
reconstructed from the bare decoupling coefficient Cas ^^ Section ^ (see Eq. (|1.5D). 



2.1 Gluon self-energy 

The bare gluon self-energy at g^ = in the full theory can be cast in the following form^ 
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TT 
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(2.1) 



where the exact dependence on e = (4 — d)/2 (d is the space-time dimension) of the bare 
two-loop result is given by 

3 + lie - £2- 15^3 _^ 4^5' 



Ph 



1 



4(2-e)(l + 2e) 



-CF-{9 + 7e-We^) + CA- 



(2.2) 



2(l-e)(3 + 2e) 

{Cp = {N^ — l)/{2Nc) and Ca = Nq are the eigenvalues of the quadratic Casimir operators 
of the fundamental and adjoint representation of SU{Nc), respectively, and Tp = 1/2 is 
the index of the fundamental representation). The three-loop quantities Phg, Phi and P}^^ 
are only available as an expansion in e. The analytic results read 
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where ^ = 1 — Cq , and fi^ 



B, = 16 Li4 ( - ) + - log2 2(log2 2-7t')- —tt' . 
\2) 3 ^ ^180 

A new result obtained in this paper is the analytic expression for Pbc{x) which arises 
from diagrams where b and c quarks are simultaneously present in the loops (see Fig. |^ for 
typical diagrams). The analytic expression is given by 



Phc{x) 



Ct. 
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Figure 1: Feynman diagrams contributing to the gluon self-energy. Thick and thin straight hnes 
correspond to b and c quarks, respectively. Wavy lines represent gluons. 




Figure 2: Feynman diagrams with two heavy-quark loops contributing to the ghost self-energy. 
The notation is adopted from Fig. |l|. 
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with L = logx, 
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where the functions L±{x) are defined in (A. 6). The function Pbc{x) satisfies the properties 



Pbcix'^) = Pbc{x) , Pfec(l) = 2Phh 



(2.5) 



which are a check of our result. For x — )• 0, the hard contribution to Pbc{x)x~^'^ is given by 
Phi- However, there is also a soft contribution, and it is not possible to obtain a relation 
between Pbdx — > 0) and Phi if they are expanded in e (this would be possible for a non-zero 

e < 0, cf. (ITgD). 



2.2 Ghost self-energy 

The bare ghost self-energy at q^ = can be cast in the form 
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where the two-loop term is given by 
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and the e expansions of the single-scale three-loop coefficients read 
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The function Cbc{x) is obtained from the diagram of Fig. |2| and can be written as 

3-2e 
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where Ilb{k'^) and nc(A;^) are the b- and c-loop contributions to the gluon self-energy. The 
integral I{x) is discussed in Appendix ^ where an analytic result is presented. In analogy 
to Eq. ( p.5D , we have 

Cbc{x-^) = Cbc{x) , Cbc{l) = 2Chh . (2.11) 



For a non-zero e < 0, Cbdx — )• 0) ^ Chix (only the hard part survives in ( |A.9| )) 



2.3 Light-quark self-energy 

The parts of the hght-quark self-energy Sy(0) and ^^(O) (with vanishing light-quark 
masses) are conveniently written in the form 
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where 
Vh 
and 



e(l + e)(3-2e) 
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Exact d-dimensional expressions for these coefficients have been obtained in |17]. 

The quantities Vbc{x) and Sbc{x) arise from diagrams similar to Fig. ^ and can be 
expressed in terms of I{x): 
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32(2 -£)''"■" """'"'' 32 

They satisfy the relations analogous to Eq. ( |2.5| ) which again serves as a welcome check of 



our calculation. Retaining only the hard part of ( A. 9 ) for x ^ 0, we reproduce Vm, Shi- 
Vbc has been calculated up to 0{e^) in Ref. [p^ ]. 

2.4 Ghost gluon vertex 

The two-loop correction vanishes in the arbitrary covariant gauge exactly in e, see Ap- 
pendix 1^. For the same reasons, the three-loop correction contains only diagrams with a 
single quark loop (bottom or charm), and vanishes in Landau gauge: 



-6e\ry2, 
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3. Decoupling for as 

The gauge parameter dependence cancels in the bare decoupling constant ( |1.5D (which 
relates a^g ^^ '^so^ ' ^^^ -^l- (L^)). Since the result is more compact we present analytical 



expressions for (Co J which reads 
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Note that Zf,c(x^ ) = Z{,c(x), ^bc(l) = 2Z/i/j. If desired, the vertices T^igq and T^iaa can 
be reconstructed using Eq. ( |1.5[ ). 

In order to relate the renormaUzed coupUngs q^ ^ (^) and as {fi), we first express ah 
bare quantities in the right-hand side of the equation 



via the MS renormahzed ones M, EO, 21 
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(and similarly for rrico). This leads to an equation where a^g is expressed via the nj- 
flavour MS renormahzed quantities^ Us ^ (/i), rndfJ') and mh{fi). In a next step we invert 
the series 
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*Note that the masses mdlJ.) and mh{^) (and mco, Jrito) are those in the full n/-flavour QCD. They 
do not exist in the low-energy nj-flavour QCD, and therefore we do not assign a superscript n/ to these 
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to express as [fi') via a^^ , and substitute the series for a^ derived above. 

In order to obtain compact formulae it is convenient to set /x = fhf, where fhh is defined 
as the root of the equation mb{fhb) = rhi,. Furthermore, we choose ^' = mc{fhb) and thus 
obtain as {mc{'mb)) as a series in as ^ {rhb) with coefficients depending on 
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are defined in such a way that dF,A{^) = 0. Thus for x = 1 Eq. ( |3.5D reduces to the ordinary 
decouphng of rih + ric flavours with the same mass Q. For x ^ 1 the functions dpix) and 
dyi(x) become 
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iTT'd^TT-b)) can be obtained by inverting the se- 
ries ( |3.5D . If one wants to express Qs"'''(/ic) as a truncated series in as ^ (/^b) (without 
resummation) for some other choice of fn, ^ rrih and fic ~ ^-c; this can be easily done in 
three steps: {i) run from fif, to m^ in the nj-flavour theory (without resummation); (ii) 
use Eq. (|3.5|) for the decoupling; and (iii) run from mc{fhb) to //c in the n;-flavour theory 



(without resummation). After that, relating as (//') and Ug ^ (/i) for any values of fi and 
fi' (possibly widely separated from m-fe and rric) can be done in a similar way: (i) run from /u 
to /if, in the nj-flavour theory (with resummation); (ii) use the decoupling relation derived 
above; and {in) run from fic to fi' in the n;-flavour theory (with resummation). The steps 
(i) and (iii) can conveniently be performed using the program RunDec [^]. 

In the case of QCD {Tp = 1/2, Ca = 3, Cp = 4/3, nb = ric = 1) the decoupling 



constant in Eq. (3^) reduces to (for e = 0) 
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For X <C 1 the coefficient of (as/vr)^ becomes 
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4. Decoupling for the light-quark masses 



The bare quark mass decoupling coefficient Cm ^^ -^1- (L^) is determined by Sy(0) and 
£5(0), see Eq. ( p. 12 ); it is gauge parameter independent. The renormahzed decouphng 
constant Cm in Eq. (|lj) (see pl| , p^ for the mass renormahzation constants) can be 
obtained by re-expressing aP in the denominator via as ^ (cf. Sect. ^; note that in Ca^ 
positive powers of e should be kept). Our result reads 



„W)/ 



.,("/)/ 



Cm(mc(mb),m6) = l+rfrgF _ ^ +^^^^^1 _ +(^'Cf\ ^ ' 



+• 



vr 



vr 



vr 



(4.1) 



where 



-^l(i-L£ + ^lV + o(.3; 



r2 r 

d- = [9Cf + HCa - ^Tpini + nc)] — - {^Cf + 97Ca - 20rF(n; + n^)] — 

o 4o 



89 ^ , 



L3 



L^ 



+ <! - \^Cf + IIC^ - ^Tpini + ne)] — + \^Cf + 97Ca - 20r^(nj + n,)\ ^^ 



+ 



37r2nb - 89n, 



72 



:.,._(..n^)T.^}..0,e^,. 



d^ 



{S^Cf + 11Ca)(9Cp + 22C^) , 27Ci. + 44Ca^ , , 
16 4 

-r|(2(n/ + nc)2 -nfeUc 



L3 



+ 



q 1S73 / 197 

-Cl + TICfCa + ^^Cl - f 9C7^ + — Ca ) TF{n, + n. 



2 20(ni + nc)2 - 29nf,nc 



9 



L2 



13 



+ 



-129CF \Cf-^^- ^^C\ - m{CF - CA)TF{ni + n,)C3 

An 1^ I')', ^^^ 11 \^^r T ^139 47 

ACfTf { 23ni + —ric - —Ub I + -C^Tp ( -^ni - —Uc - Suf, 



27 



r| ((n; + nc)i35ni + 124nc) + Uirihiic 



L 

64 



+ 



+ 



X V^ " 20 + ¥^=^ " 144 j " ^ V 10 ^ ^^' 



16627\ 

1944 y 



-5 



2\2 



(1-X^) 



■(l2 + 2) -2 



1-x^ 



L 



96 



+ 0(e). 



At a; = 1 this result reduces to the ordinary decoupling of n^ + ric flavours with the same 
mass 10. 

Specifying to QCD leads to (for e = 0) 
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5. Decoupling for the fields 

5.1 Gluon field and the gauge parameter 

Decoupling of the gluon field and the gauge fixing parameter are given by the same quantity 

(5.1) 






In a first step we replace the bare quantities in the right-hand side via the renormalized 
ones using Eqs. (O). (Fl)> and p|, |2l, |0| 
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and thus we express Oq via the nj-flavour renormalized quantities. In a next step we can 

(5.3) 



find a'"'^(/i') in terms of Oq by solving the equation 
iteratively. The result reads 
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with a = a^^f>{mb). The easiest way to express a^"'f'{mb) via a^'^^\mc{mb)) is to re- 
express Us ■' (m-b) via as {vac{fhb)) in the right-hand side of the equation a^"-'' {mc{mb)) = 
a^"'f\mb)CA{'ffib,'>TT,c{'mb)) and then solve it for a^"'f'{fhb) iteratively. 



5.2 Light-quark fields 

The bare decoupling coefficient Cq of Eq. (|1.2| ) is determined by Sy (0) (cf. Eq. ( p.l2D ). The 
renormalized version (q (|1.7D can be obtained (see Refs. [24, ^, ^ for the three- loop wave 
function renormalization constant) by re-expressing Os and a^"'^' in the denominator via 
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the nj-flavour quantities (see Sects. ^ and 5A; note that positive powers of e should be 
kept). The result can be cast in the form 
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Note that the power corrections in x drop out in the sum of all diagrams. For x = 1 this 
result reduces to the ordinary decoupling of n^ + Uc flavours with the same mass |0] (see 



Ref. [17] for an expression in terms of Ca and Cf)- 

5.3 Ghost field 

The bare decoupling coefficient Cc ™ ^1- (|1-2| ) is determined by llc{0) as given in Eq. (p.6|). 
The renormalized decoupling constant Cc of Eq. ( |1.7| ) is given by (see Refs. |25, 20 1 for the 
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corresponding renormalization constant) 
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Figure 3: ct-s (Mz) as obtained from ots {Mr) as a function ^dcc. The dashed lines (long 
dashes include higher order perturbative results) correspond to the single-step approach and the 
dash-dotted curves (short dashes: /idcc,c — A*dcc, long dashes: ^dcc,b — /^dec) are obtained in the 
conventional analysis using four-loop running and three-loop decoupling relations. The dotted line 
results from a five-loop analysis of the two-step (see text for details). 



6. Phenomenological applications 

In this section we study the numerical consequences of the decoupling relations computed 
in the previous sections. For convenience we use in this Section the decouphng relations 
in terms of on-shell heavy quark masses (see Appendix ^ and the Mathematica file which 
can be downloaded from [E^]) which we denote by Mc and M^. 



6.1 JI'\Mz) 



from aF(M^) 



Let us in a first step check the dependence on the decouphng scales which should become 
weaker after including higher order perturbative corrections. We consider the relation bet- 
ween Us (M-r) and as (Mz)- cis (Afr) has been extracted from experimental data using 
perturbative results up to order a^ p^. Thus it is mandatory to perform the transition 
from the low to the high scale with the highest possible precision. In the following we 
compare the conventional approach with the single-step decoupling up to three-loop order. 

For our analysis we use for convenience the decoupling constants expressed in terms of 
on-shell quark masses. In this way the mass values are fixed and they are not affected by 
the running from Mr to Mz- In our analysis we use Mc = 1.65 GeV and Mb = 4.7 GeV. 
Furthermore, ai (M,-) = 0.332 ^^ is used as starting value of our analysis. 

In Fig. |^(a) we show ai (Mz) as a function of /idee, the scale where the c and b quarks 



(3), 



„(3), 



are simultaneously integrated out. In a first step a^ (M,-) is evolved to a^ (^dec) using the 
A^-loop renormalization group equations. Afterwards the {N — l)-loop decoupling relation 
is applied and finally A^-loop running is employed in order to arrive at al (Mz)- One 
observes a strong dependence on //dec for A^ = 1 (short-dashed line) which becomes rapidly 
weaker when increasing N leading to a reasonably flat curve for A = 4 (longer dashes 
correspond to larger values of A). 
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6.2 Comparison of one- and two-step decoupling approach 

In the step-by-step decoupling approach we have two decoupling scales fJ-dec,c and fJ-decfi 
which can be chosen independently. First we choose^ /^dec,c = 3 GeV and identify /Udec.fe 
with /idee- The result for A^ = 4 is shown in Fig. ^b) together with the four- loop curve from 
Fig. ^(a) as dash-dotted line (long dashes). One observes a significantly flatter behaviour 
as for the one-step decoupling which can be explained by the occurrence of log(^^/M^) 
terms in the one-step formula which might become large for large values of /i = /idee- 
Alternatively it is also possible to study the dependence on fj.dec,c, i-e., identify /idec,c with 
/Xdec) set /idee, 6 = 10 GcV and compare to the one-step decoupling. The results are also 
shown in Fig. |^(b) as dash-dotted line (short dashes) where only values /idee ^ 10 GeV are 
considered. 

For comparison we show in Fig. y(b) also the result of the two-step five-loop analysis 
as dotted line where the four-loop decoupling relation is taken from Refs. ^, |5|. The 
(unknown) five-loop coefficient of the f3 function, /34, is set to zero.^ If one restricts to 
scales /fdec between 2 GeV and 10 GeV it seems that the four-loop decoupling constant 
is numerically more relevant than the power-suppressed terms included by construction in 
the one-step decoupling procedure. Thus, from these considerations one tends to prefer 
the two-step decoupling over the one-step approach as it seems that the resummation of 
log(/i^/M^j,) is more important than the inclusion of power-suppressed corrections. 

Let us in a next step restrict ourselves to decoupling scales which are of the order of 
the respective quark masses. In Tab. Ill we compare the value for a^ (Mz) as obtained 
from the one- and two-step decoupling where two variants of the former are used: Cos 
which directly relates Us {fJ-c) and ai (/*&) as given in Eq. (1.7) with jj' = fic and fi = fi/, 



(C«s(Mci/ife); see also ||27|]) and the version with only one decoupling scale where /i' = /i has 
been set (^^^(/i)). We thus define two deviations 

6a'f^ = afHMz) ^ ^ - af\Mz) 

Cas(t^c,M 2-step 

5ai^) = af)(Mz) -af)(Mz) , (6.1) 



Cas (m) 



2-step 



where the scale fj, in the second equation is either identified with fic (right part of Tab. |l]) 
or /if, (left part), respectively. 

It is interesting to note that (except for the choice /tc = 2 GeV and //;, = 10 GeV) 
the deviations presented in Tab. || amount to about 30% to 50% of the uncertainty of the 



world average for as{Mz) which is given by 6as = 0.7 • 10 [30| 



6.3 Improving the t^vo-step approach by power-suppressed terms 

From the previous considerations it is evident that the resummation of logarithms of the 
form [as log(/ic//ife)] , which is automatically incorporated in the two-step approach, is nu- 
merically more important than power-suppressed terms in Mc/Mb. Thus it is natural to use 



^It has been argumented in Refs. PSJ] that in the case of charm the scale fj, = ttIc is too small leading to 
a value of as which is too large. Thus mc(3 GeV) has been proposed as reference value. 
®For /34 > the dotted curve in Fig. Mb) moves towards the four-loop curve. 
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I^b 


«F(Mz) 


r (a) 


6af^ 


(GeV) 




xlO^ 


xlO=^ 


2 


0.11985 


-0.28 


0.18 


5 


0.11977 


0.23 


-0.16 


7 


0.11974 


0.36 


-0.26 


10 


0.11970 


0.19 


-0.36 



fJ-c 


af\Mz) 


X (a) 

Sal 


Sa^^ 


(GeV) 




xl03 


Xl03 

(M = /^c) 


2 


0.11984 


-4.02 


0.20 


3 


0.11970 


0.19 


0.14 


4 


0.11961 


0.33 


0.10 


5 


0.11955 


0.26 


0.06 



(5) 

Table 1: Decoupling scale ah (Mz) as obtained from the four- loop analysis of the two-step 
approach, and the deviations as defined in the text. In the left table /ic = 3 GeV and in the right 
one /Zb = 10 GeV has been chosen. 



the two-step approach as default method and add the power-corrections afterwards. This 
is achieved in the following way: In a first step we invert Casil^cfJ-b) (cf- Eq. (|1.7D) and ex- 



press it in terms of as (/Uc) in order to arrive at the equation as (idb) = Ca/ (/"c l^b)as (/i, 
Now an expansion is performed for Mc/Mf, — )• to obtain the leading term which is then 
subtracted from Cq/(/"c> /"{,) since it is part of the two-step decoupling procedure. The 
result is independent of fic and fib and has following series expansion 



5c; 




6661 1409 , 1 .V , V ,^, ,, 
18000 21600 160 ' ^ ' 



0.170 



(6.2) 



where the numerical value in the second line has been obtained with the help of the exact 
dependence on x. Note that the linear term in x arises from the MS-on-shell quark mass 
relation. The quantity 5C,a^ is used in order to compute an additional contribution to 
al ifJ-b) as obtained from the two-step method: 



fef(^,) = <5C->f)(^^) 



(6.3) 



Inserting numerical values leads to shifts which are at most a few times 10~^ and are thus 
beyond the current level of accuracy. It is in particular more than an order of magnitude 
smaller than the four- loop decoupling term which is shown as dotted curve in Fig. ^(b). 

Note that as far as the strong coupling in Eq. ( |6.2D is concerned both the number of 
flavours and the renormalization scale of q^ are not fixed since power-suppressed terms 
appear for the first time at this order. However, the smallness of the contribution is not 
affected by the choices made in Eq. (6.2). 



6.4 One-step decoupling of the bottom quark with finite charm quark mass 

An alternative approach to implement power-suppressed corrections in rric/mb in the de- 
coupling procedure is as follows: We consider the step-by-step decoupling and use at the 
scale //dec c the standard formalism for the decoupling of the charm quark as implemented 
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in RunDec ||23[|. At the scale /idec.b) however, we consider the matching of five- to four- 
flavour QCD where we keep the charm quark massive. This requires a modification of the 
formulae in Eqs. ( |1.2| ) and ( |1.5| ) to {n'r = rij — 1) 

where the nj-flavour quantities contain contributions form massive charm and bottom 
quarks. They are identical to the one-step decoupling procedure described above. In the 
n'r-flavour quantities appearing in the denominators those diagrams have to be considered 
which contain a charm quark. Note that they depend on the bare parameters of the 
effective theory (a^Q^ , Oq ^ , m^Q^ ) and thus they have to be decoupled iteratively in 
order to express all quantities on the r.h.s. of the above equations by the same parameters 
(a^Q-^ , Oq ^ , m^^ ). In the standard approach the n'/-flavour quantities vanish since only 
scale- less integrals are involved. 

As a cross check we have verified that we reobtain the analytical result for the single- 
step decoupling if we apply the formalism of Eq. ( |6.4D and the subsequent decoupling of 
the charm quark at the same scale. 

We have incorporated the finite charm quark mass effects in the two-step decoupling 
approach (cf. Fig. ^ and observe small numerical effects. A minor deviation from the 
nic = curve can only be seen for decoupling scales of the order of 1 GeV which confirms 
the conclusions reached above that the power-suppressed terms are numerically negligible. 
Thus we both refrain from explicitly presenting numerical results and analytical formulae 
for the renormalized decoupling coefficients as obtained from Eqs. (|6.4D. 



6.5 Decoupling effects in the strange quark mass 

In analogy to the strong coupling we study in the following the relation of the strange 
quark mass ^^(/x) defined with three and five active quark flavours, respectively. The 
numerical analysis follows closely the one for a^: A^-loop running is accompanied by (N — 
l)-loop decoupling relations. It is, however, slightly more involved since besides ms{fj.) 
also as(/i) has to be known for the respective renormalization scale and number of active 
flavours. We organized the calculation in such a way that we simultaneously solve the 
renormalization group equations for ms{fi) and as(//) (truncated to the considered order) 
using Mathematica. 

In Fig. H we show rris {Mz) as a function of /^doc and again compare the single-step 
(dashed lines) to the two-step (dash-dotted lines) approach. For our numerical analysis 
we use in addition to the parameters specifled above ms{2 GeV) = 100 MeV. The same 
conclusion as for ag can be drawn: The difference between the two approaches becomes 
smaller with increasing loop order. At the same time the prediction for rris {Mz) becomes 
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N 



in 



0.064 



0.062 = 



0.06 



!/3 



0.058 



0.056 




1 



10 



Figure 4: 



,(5) 



(Mz) as a function of /idee- The dashed lines correspond to the single-step approach 



and the dash-dotted curves are obtained in the conventional analysis (with /idoc,c — 3 GeV and 



A'dc 



/idee)- Longer dashes correspond to higher loop orders. See text for more details. 



more and more independent of ^doc- The results again suggest that the power-corrections 
Mc/Mf, are small justifying the application of the two-step decoupling. 

7. Effective coupling of the Higgs boson to gluons 

The production and decay of an intermediate-mass Higgs boson can be described to good 
accuracy by an effective Lagrange density where the top quark is integrated out. It contains 
an effective coupling of the Higgs boson to gluons given by 



-^cff 



-CiOi 



(7.1) 



with Oi = G^yG^^ . Ci is the coefficient function containing the remnant contributions of 
the top quark, G^^'^ is the gluon field strength tensor, cj) denotes the CP-even Higgs boson 
field and v is the vacuum expectation value. 

The effective Lagrange density in Eq. ( [7.1| ) can also be used for theories beyond the 
Standard Model like supersymmetric models or extensions with further generations of heavy 
quarks. In all cases the effect of the heavy particles is contained in the coefficient function 
Ci. 
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In Ref. [^ a low-energy theorem has been derived which relates the effective Higgs- 
gluon coupling Ci to the decoupling constant for Og ■ In this Section we apply this theorem 
to an extension of the Standard Model containing additional heavy quarks which couple 
to the Higgs boson via a top quark-like Yukawa coupling. Restating Eq. (39) of Ref. [^ in 
our notation and for the case of several heavy quarks leads to 



C, 



-I Nh 



i=l 



dM. 



log Co 



(7.2) 



where Nh is the number of heavy quarks with on-shell masses Mi . Using Ca^ from Eq. ( |1.7| ) 
(see also p^) we obtain for Ci the following result'^ 



Ci 




«.r^(/^) 



5 
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CaTf^ + T'^-t^ 



18 



Nh 



+ 



' (full)/ \' 

Oil (/i) 



vr 



2 9 
-CpTp-—Nh + CfCaTf 
d4 



72 96 



+ CfTI 
+CaTI 



5 17 o 

^^N,n, + —Nl 






CITf 



7 



1063^^ 



47 49 ^, 5 ^ ■ 

n; iV;, H l^h 

864 1728 24 



Nh 



rp3 v2 ^^ 



'F^ft 



54 



(7.3) 



(full) 



where al is the strong coupling in the full theory with ni + Nh active quark flavours 
andS, = E£'ilog(/xVM^) 



After expressing a^ in terms of al and specifying the 



colour factors to SU(3) we reproduce the result of Ref. |31| which has been obtained by an 
explicit calculation of the Higgs-gluon vertex corrections. For A'^^ = 1 the result obtained 
in Ref. [^J is reproduced. It is remarkable that although Ca^ contains di- and tri-logarithms 
there are only linear logarithms present in Ci. 



8. Conclusion 

The main result of this paper is the computation of a decoupling constant relating the 
strong coupling defined with three active flavours to the one in the five-flavour theory. At 
three-loop order Feynman diagrams with two mass scales, the charm and the bottom quark 
mass, have to be considered. The corresponding integrals have been evaluated exactly and 
analytical results have been presented. The new results can be used in order to study 
the effect of power-suppressed terms in Mc/Mh which are neglected in the conventional 
approach j^. Various analyses are performed which indicate that the mass corrections 
present in the one-step approach are small as compared to log(/i^/M^^) which are resummed 
using the conventional two-step procedure. 

Using a well-known low-energy theorem H we can use our result for the decoupling 
constant in order to obtain the effective gluon-Higgs boson coupling for models containing 



''Note that up to three-loop order there are only diagrams with at most two different quark flavours. 
Thus it is possible to obtain the result for Ci for Nh heavy quarks. 
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Figure 5: Master integral /i with four massive lines. Thick and thin straight lines correspond to 
b and c quarks, respectively. Master integral I2 contains an additional numerator. 



several heavy quarks which couple to the Higgs boson via the same mechanism as the top 
quark. This constitutes a first independent check of the result presented in Ref. [^] where 
the matching coefficient has been obtained by a direct evaluation of the Higgs-gluon-gluon 
vertex diagrams. 
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A. Integral I{x) 



With the help of FIRE [13| we can express the integral I{x) as defined in Eq. ( 2.1[1| ) as a 



linear combination of master integrals 

lix) = /(x^i) = (^_i)(^_4)(^_6)(^_8)(^_io) ^^-^^ 

X - (cio + cu{x-^ + x^) + cuix'^ + x^)) hix) 

+ — (d - 2)(X-1 + X) (C20 + C2i{x-^ + x^)) hix) 
lb 

c_i(x2+^ + a;"2-E) + co(x^ + x"^) + ci(x"^+^ + x^'") + C2(x~''+^ +x*^~^)' 
(d-2)2(d-3)((i-5)(d-7) 

/i and I2 are master integrals with four massive lines (see Fig. ^) which are given by 



/i(x)=/i(x-i) 
h(x) = hix 



(mfemc)-2+3^ /■ d'^ki d'^k2 d'^h 



(i^d/2)3r3(e) J DiD2D^Di ' 
1, _ {mbmc)-^^^^ f N d'^ki d'^k2 d'^k^ 



(m««/2)3r3(e) J DiD2DiDi 
h = 'ml-kl, D-i = ml-h 
Di = ml-{ki-k2 + k^f , N=-{ki-k2f, (A.2) 



Di = ml - kj , D2 = ml-kl, D3 = ml - kj , 
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and Cj and Cij are coefficients depending on d = 4 — 2e 

cio = (d - l)(5(i^ - 104(i3 + 73^^ - 2116d + 2086) , 

cii = {d- l)(2(i - 7)(2(i^ - 35^^ + 180d - 256) , 

ci2 = (d - 9)(2(i - 5){2d - 7){2d - 9) , 

C20 = 2((i^ - 22d^ + 165^2 - 491d + 487) , 

C2i = (rf-9)(2(i-7)(2(i-9), 

c_i = (d - 3)(d - 5)(d - 7)((i - 9)(2d - b){2d - 7){2d - 9) , 

co = {d- l){d - 3){id^ - W8d^ + 1090d^ - 5009^^ + 9838d - 5335) 

ci = {d- l){d - 7){2d^ - 46d^ + 384d^ - 1423^2 + 2158d - 739) , 

C2 = {d- l){d -5){d- 7){d - 9){2d - 7){2d - 9) . 

Tlie master integrals used in Ref. [fl] are related to /i 2 by 



/4.3 = {mbmcy'''T'{e)Ii{x), 



h.3a = {mbm,y-''r\e) 



1 — x^ 



-\{d-3- {2d - 5)x^) h{x) + ^{d - 2)xh{x) + — 



x^ + x'^ ^ 



Using their expansions in e 113] we obtain 






^. + i(f + 3L^).^ + i.e^ + 



where 



32 „ _ ,, (l + x2)(5 + 22x2 + 5x^) , ,, 

-B = 64L+ X - ^ '-^ 3 '-L.{x) 

3 X'^ 



+ 



5 + 18x2 + 5x4 



L' 



10- 



1-x^ 



-L + 10- 



fl-x 



2\2 



X^ 



X^ 



64^ 

+ yC3 



1256 
81 



and 



L±{x) = L±{x ^ 
± 



Li3(x) -LLi2(x)- Ylog(l-x) + — 



^2 J?. 

Li3(-x) - L Li2(-x) - — log(l +x) + — 



with L = logx. Note that the functions L±{x) are analytical from to +cx3. 
For X = 1, 12(1) is not independent |T^ : 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



The expansion of /i(l) in e has been studied in Refs. Il^, p^]. 



(A.7) 
Using the explicit formu- 



las (3.2) and (2.3) from [14|, it is easy to get 



32 

'27 



2 25 o 

1 - -eH e^ 

3 6 



rc, + ^|e^ + 



(A.^ 
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in agreement with 

For a; — )• 0, two regions [^ contribute to I{x) (see Eq. ( 2.10[ )), the hard (A; ~ mf,) and 
and the soft {k ~ rric) one. The result for the leading term is given by 

I{x) = hx^' [1 + 0{x^)] + Isx-' [1 + 0{x^)] , (A.9) 

d-5 r(l-e)r2(l + 2e)r(l + 3e) 






3 (d - l){d - 3){2d - 9){2d - 11) 
d-6 



3{d-2){d-5){d-7) ■ 
Expanding this formula in e we reproduce Eq. (| 



r2(l + e)r(l + 4e) 



for X — )• 0. 



B. Ghost— gluon vertex at two loops 

We need this vertex expanded in the external momenta up to the linear terms. Let us 
consider the right-most vertex on the ghost line: 



A^^^Pu. 



The tensor A^'^ may be calculated at zero external momenta, hence A^'^ = Ag^^ . Therefore 
all loop diagrams have the Lorentz structure of the tree vertex, as expected. 
Now let us consider the left-most vertex: 




k.M 



M^- 



It gives k \ thus singling out the longitudinal part of the gluon propagator. Therefore, all 
loop corrections vanish in Landau gauge. Furthermore, diagrams with self-energy insertions 
into the left-most gluon propagator vanish in any covariant gauge: 



>-*-:^ 





In the diagrams including a quark triangle, the contraction of k^ transfers the gluon 
propagator to a spin propagator and a factor k^ which contracts the quark-gluon vertex. 
After decomposing ^ into a difference of the involved fermion denominators one obtains 
in graphical form 



ao 
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ao 




- -/-> J- - 




- ^- -> - -t»- - 



The diagrams with a massless triangle vanish. The non-vanishing diagrams contain the 
same Feynman integral, but differ by the order of the colour matrices along the quark line, 
thus leading to a commutator of two Gell-Mann matrices. 

The remaining diagram contains a three-gluon vertex with a self energy inserted in the 
right-most gluon propagator. The contraction of k^ with the three-gluon vertex cancels 
the gluon propagator to the right of the three-gluon vertex: 





-> ^- - - ^- -> 

The colour structure of the three-gluon vertex is identical to the commutator above, how- 
ever with opposite sign. Therefore, after summing all contributions the result is zero. 



C. Decoupling at on-shell masses 

For some applications it is convenient to parametrize the decoupling constants in terms of 
the on-shell instead of MS quark masses. The corresponding counterterm relation reads 



rribo 



os{nf) { (rif) 






a 



sQ 



Mb 



ruco 



os{ny) 



Q 



sO 



Mr 



(C.l) 



where in our application Z'., 
been calculated in Ref. 

■os{nf) 



os(nf) 



nil, 



and Z, 



os(nj) 



are needed to two-loop accuracy. They have 



zosirif) 
'mi, 

in Ref. 



(see also 



12|). Note that the two- loop coefficients of 



Znii, and z, 



are non-trivial functions of mc/mb; a compact expression can be found 



m 



The advantage of using on-shell masses is that they are identical in all theories (with 
any number of flavours). Furthermore their numerical value does not depend on the renor- 
malization scale. However, it is well known that usually the coefficients of perturbative 
series for physical quantities grow fast when expressed via on-shell quark masses and hence 
the ambiguities of the mass values (extracted from those observable quantities) are quite 
large. Nevertheless, using on-shell masses in intermediate theoretical formulae (at any 
finite order of perturbation theory) can be convenient. 

The decoupling relations are particularly compact if as (Mc) is expressed as a series 
in tts ^ (Mb) since then the coefficients only depend on Xqs = Mc/Mb (see results in [p^]). 
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